The TeV spectrum of H 1426+428 
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ABSTRACT 

The BL Lac object H1426+428 was recently detected as a high energy 7-ray source by the VERITAS 
collaboration (Horan et al. 2002). We have reanalyzed the 2001 portion of the data used in the detection 
in order to examine the spectrum of H1426+428 above 250 GeV. We find that the time-averaged spectrum 
agrees with a power law of the shape 
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The statistical evidence from our data for emission above 2.5 TeV is 2.6 a. With 95% c.L, the integral 
flux of H1426+428 above 2.5 TeV is larger than 3% of the corresponding flux from the Crab Nebula. The 
spectrum is consistent with the (non-contemporaneous) measurement by Aharonian et al. (2002) both in 
shape and in normalization. Below 800 GeV, the data clearly favours a spectrum steeper than that of any 
other TeV Blazar observed so far indicating a difference in the processes involved either at the source or 
in the intervening space. 

Subject headings: BL Lacertae objects: individual (H 1426+428) — 7 rays: observations 
1. Introduction 

The BL Lac object H1426+428 was discovered at 
optical wavelengths with a redshift of 0.129 by Remil- 
lard et al. (1989). Measurements of the spectral en- 
ergy distribution of this close object are still sparse and 
do not sufficiently cover the expected broad two-peak 
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structure that is familiar from other BL Lac objects 
(see, e.g., Donato et al. (2001)). The first peak is ex- 
pected at X-ray energies and is thought to represent 
the synchrotron emission from relativistic electrons in 
the source. The second peak is expected at 7-ray en- 
ergies and is explained in so-called leptonic models as 
stemming from low-energy photons which are inverse- 
Compton scattered to 7-ray energies by the same pop- 
ulation of relativistic electrons that causes the syn- 
chrotron radiation. Alternative models attribute the 
second peak and also part of the X-ray emission to pro- 
cesses involving protons which are co-accelerated with 
the X-ray-emitting electrons or even partially produce 
them. For recent reviews see, e.g., Ghisellini (2000), 
Rachen (2000), Sikora & Madejski (2000). 

A lower limit on the position of the X-ray peak of 
H1426+428 was placed by Costamante et al. (2001) at 
100 keV. Currently the only published measurements 
above 100 keV come from Cherenkov telescopes above 
280 GeV (Horan et al. 2002) and above 700 GeV (Aha- 
ronian ct al. 2002). Costamante & Ghisellini (2002) 
predict the peak of the 7-ray emission of H1426+428 
to be at several ten GeV. 

The study of the high-energy spectrum of this "new" 
TeV source is especially interesting since its rcdshift is 
four times as large as that of Mkn 421 and Mkn 501, the 
only other extragalactic objects detected at TeV ener- 
gies with good spectral information. Due to this larger 
distance, it is expected that signs of absorption of the 
7-radiation via interaction with the intergalactic opti- 
cal and infra-red background (IIRB) will be more pro- 
nounced, possibly permitting one to infer constraints 
on the IIRB photon density. 

In this letter, we present the results of a spectral 
reanalysis of H1426+428 observations made with the 
Whipple 10 meter 7-ray telescope on Mt. Hopkins, 
Arizona, in the first half of 2001. We use our standard 
method described in Mohanty et al. (1998) for deriving 
the spectrum. Due to the weakness of the source and 
the special interest in the emission at the highest ener- 
gies, we then examine the emission above 2.5 TeV with 
specially developed cuts. 

2. Dataset 

The detection of H1426+428 (H1426 in the follow- 
ing) as presented in Horan et al. (2002) is based on 
four separate datasets from 1995-98, 1999, 2000 and 
2001 respectively. Since the sensitivity of the Whipple 
telescope was considerably improved over from 1995 to 



Table 1: General properties of the H1426 dataset and 
the corresponding OFF dataset. 





ON 


OFF 


first MJD 


51940 


51875 


last MJD 


52275 


52275 


raw number of events 


3698176 


3826088 


min. zenith angle 


10° 


6° 


max. zenith angle 


38° 


38° 


mean zenith angle 


18.0° 


18.7° 


total observation time 


38.1 h 


39.9 h 



2001 and the largest fraction of the observation time 
was invested in 2001, it was only in the 2001 dataset 
that a 7-ray signal with significance > 5tr was found. 
For this observing period, contemporaneous observa- 
tions of the Crab Nebula are available (Krennrich et al. 
2001) containing high-significance 7-ray signals which 
can serve to calibrate the telescope. In order to de- 
rive a spectrum for H1426, we have therefore chosen to 
confine ourselves to the 2001 dataset. 

The dataset used here consists of 87 pairs of 
on-source (ON) and off-source (OFF) observations 
("runs"). The run duration for the ON data varied 
between 10 and 28 minutes resulting in a total expo- 
sure time of 38.1 hours. The OFF dataset serves to 
determine the background caused by hadronic cosmic 
rays. Out of the 87 ON-OFF pairs in our dataset, 34 
were taken in actual ON-OFF mode, i.e. with exactly 
matching zenith angle distribution and night sky back- 
ground conditions. The OFF data for the remaining 46 
ON runs was obtained by selecting matching OFF runs 
out of the pool of other OFF data taken in early 2001. 
The matching was based on how close to each other the 
runs were in terms of date, weather conditions, eleva- 
tion and night sky brightness. This procedure results 
in a good agreement of the zenith angle distribution, 
but small discrepancies remain. Also, some of the OFF 
runs have a longer duration than their ON partners. As 
a result, the total OFF exposure time is 4.7 % longer 
than the ON exposure time. These differences were 
eliminated in later stages of the analysis by using the 
number of events not coming from the source direction 
to scale the background rate. This scaled background 
should be equal for both the ON and the OFF data. 
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3. Data analysis 

3.1. Determination of image parameters 

The determination of the image parameters WIDTH, 
LENGTH, ALPHA, SIZE etc. for ON and OFF data 
was carried out following the standard procedure de- 
scribed in Reynolds et al. (1993) and Cawley (1993). 

3.2. Monte Carlo simulation 

Our data was taken between 10° and 38° zenith an- 
gle with less than 10% taken at zenith angles between 
30° and 38°. The average zenith angle is near 20°. The 
properties of air showers change as a function of the 
cosine of the zenith angle. Given the relatively large 
statistical error of the gamma signal from H1426 (see 
below) , adequate accuracy for the determination of the 
detector sensitivity can therefore already be reached by 
simulating air showers only at the average zenith angle. 

The shower simulation code ISUSIM (based on "kas- 
cade" by Kertzman & Sembroski (1994), see Mohanty 
et al. (1998)) was used to generate 3.8 x 10 6 show- 
ers at zenith angle 20° between 0.05 TeV and 100 TcV 
with a spectral index of 2.5 (MC set 1). Other spectral 
indices were simulated by re-weighting. The simula- 
tion of the Whipple telescope was carried out with the 
"simveritas" code (Duke & LcBohec 2002, VERITAS 
internal report) . The telescope parameters in this sim- 
ulation were identical with those used in Krennrich et 
al. (2001). 

3.3. Derivation of the "Extended Cuts" 

Using MC set 1 and the spectrum analysis algo- 
rithms developed in Mohanty et al. (1998), the "ex- 
tended cuts" for the 7-hadron separation in our dataset 
were derived. These utilize only the image parame- 
ters WIDTH, LENGTH and ALPHA. They are derived 
as a second-order polynomial function of log(SIZE) 
by approximating the distributions of WIDTH and 
LENGTH and ALPHA as Gaussians and determining 
the mean a and the standard deviation a as a function 
of log(SIZE). The lower cuts on WIDTH and LENGTH 
respectively are then of the form (a — ta) while the up- 
per cuts are of the form (a + ta). The upper cut on 
ALPHA is also of the latter form. The constant t, the 
"tolerance", determines the efficiency of the cuts. It 
is generally chosen to be 2.0 in order to minimize pos- 
sible systematic effects arrising from imperfections of 
the Monte Carlo simulations. In addition to the ex- 
tended cuts, standard quality cuts were imposed as in 



Krennrich et al. (2001). 

3.4. Background determination 

Due to the 5 % difference in exposure time and the 
remaining differences in elevation and night sky back- 
ground, the number of events in the region ALPHA> 
30° (where according to simulations there should be 
no significant contribution caused by 7-rays from the 
source) is smaller for the ON data (184440 events) than 
it is for the OFF data (203531 events). The ratio r of 
these numbers (here 0.9062) can be used to scale the 
number of OFF events at small ALPHA in order to 
determine the background. Figure 2a shows the scaled 
OFF data distribution (with r = 0.9062) superimposed 
on the ON distribution. 

After applying the extended ALPHA cut (see Fig- 
ure 1) the number of events is 51819 ON and 50279 
OFF with an excess of 1540 events and significance 5 
a. Figure 2b shows that the shape of the ON and the 
OFF distributions agrees well in the normalization re- 
gion (ALPHA> 30°). 

3.5. Energy reconstruction and binning 

An algorithm described in Mohanty et al. (1998) was 
used to derive the coefficients in an analytical expres- 
sion for the estimated primary 7-ray energy E est from 
the Monte Carlo dataset 1 (see above). This expression 
is of the form 

log 10 (-E es t) = ai + a 2 L + a 3 D + a 4 L 2 + a 5 D 2 + a 6 DL 

where L = log 10 (SIZE), D = DIST. This estimator 
gives a roughly constant energy resolution above 300 
GeV. The energy resolution (the standard deviation in 
log 10 (i?ost)) was determined near 400 GeV to be 0.166, 
i.e. AE/E = 33 % . This value is used to determine the 
binning of the data with respect to log 10 (-E es t)- Statis- 
tics theory (Scott 1979) suggests that the ideal compro- 
mise between maximal structure resolution and mini- 
mal error in spectral parameters is achieved for weak 
signals if the energy bin width is chosen to be two stan- 
dard deviations. In our case we must choose the energy 
bin width equal to 2 x 0.166 = 0.332. 

The collection area of the telescope after our cuts 
on the image parameters reaches 10 % of the maximum 
value (75,000 m 2 ) at an energy of 250 GeV. The peak 
7-rate for a standard spectrum of index 2.7 is observed 
near 410 GeV. We choose 250 GeV to be the lower 
edge of the second energy bin. The first energy bin 
will be excluded from fits of spectral functions. We 
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Fig. 1. — The upper (solid) and lower (dashed) 
"extended cuts" on the image parameters WIDTH, 
LENGTH and ALPHA as a function of log(SIZE) with 
tolerance t — 2.0. These cuts minimize systematic er- 
rors stemming from the simulations. The dotted lines 
show the upper cuts for maximum significance above 
2.5 TeV (optimized on Crab Nebula data, see section 
3.7). 



find this necessary since the systematic errors of the 
simulation at the trigger threshold are much larger than 
at higher energies. The first column of table 2 gives the 
logarithmic centers of the equidistant energy bins which 
result from these choices. 

In order to ensure that the normalization constant 
r is independent of E est , the ratio was measured sepa- 
rately in each energy bin. We find that the estimated 
energies in the H1426 ON dataset are systematically 
smaller than those in the OFF dataset by a small con- 
stant factor. By varying the factor until the scaling 
ratio r is the same within the statistical errors in all 
energy bins, we determine the value of K, the neces- 
sary correction factor for the estimated energies in the 
OFF data, to be K = 0.9784 ± 0.005. 

This small effect and its correction were studied in 
detail with datasets from observations of the Crab Neb- 
ula and Mkn 421 (two strong sources) and was found to 
be a result of the differences in the nightsky background 
conditions and zenith angle distributions. It was veri- 
fied that the application of the correction factor K to 
the estimated energies of the OFF dataset results in 
correct spectral measurements. The statistical error of 
K leads to an increase of the statistical errors of the 
spectral parameters as discussed below. 
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Fig. 2.— (a) (left) ALPHA distribution for the H1426 
dataset ON (solid) and OFF (dashed, scaled by a fac- 
tor 0.9062) after application of the extended cuts with 
t = 2.0 described in section 3.3 except the cut on AL- 
PHA, (b) (right) the difference between the two distri- 
butions shown in (a) with a constant function fitted to 
the region ALPHA> 30°. 



3.6. Derivation of the H1426 spectrum 

After application of the correction to the estimated 
energy and the scaling to the OFF data, we bin both 
ON and OFF events after all cuts in log 10 (-E eBt ). For 
each of the energy bins, we determine a modified col- 
lection area A a . Its value depends on the assumed 
spectral index a since the collection area has to be av- 
eraged over a finite energy bin width and the effects 
of the finite energy resolution will lead to spill-overs 
mostly from lower to higher bins, especially if the spec- 
trum is steep. The differential flux in the energy bin i 
is then determined by 



(-) 



ON, - OFF; 
A~TAEi 



■ k 



where T is the observation time and AEi is the energy 
bin width. The factor k corrects for the approximate 
calculation of the energy derivative over a finite bin 
width. It has for our case the value 0.9313 . Table 2 
shows the event statistics for each energy bin. 

Starting with the assumption that the spectral in- 
dex will be similar to 2.5, we calculate the differential 
fluxes in each energy bin using A 2 .5 and fit a power 
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Table 2: The event statistics for the 2001 H1426+428 
dataset after application of all extended cuts (tolerance 
t = 2.0). 



Energy 


ON 


OFF 


ON-OFF 


S 


(TeV) 


(events) 


(events) 


(events) 


(') 


0.171 


10341 


10454 


-113 ± 144 


-0.78 


0.369 


23149 


22255 


894 ± 213 


4.2 


0.794 


13887 


13294 


593 ± 164 


3.6 


1.71 


3424 


3283 


141 ± 81 


1.7 


3.69 


776 


760 


16 ± 39 


0.41 


7.94 


209 


200 


9 ± 20 


0.46 


17.1 


33 


33 


± 8 


0.0 



law to the resulting values excluding the bin 1 with its 
large systematic errors. After iterating twice, using ev- 
ery time the spectral index a resulting from the fit to 
recalculate A a and then recalculating the differential 
fluxes and repeating the fit, we reach convergence at 
a = 3.50 ±0.30. 

If we vary the correction K of the estimated energy 
for the OFF data within its statistical errors, we find 
that the resulting a varies by ±0.18. This has to be 
regarded as an additional, to a good approximation un- 
corrected, error since the uncertainty of K stems from 
a finite number of events which are different from those 
events used to measure a. The total statistical error on 
a is therefore V0.30 2 + 0.18 2 = 0.35. The systematic 
error on a for our method has to be determined from a 
dataset with small statistical errors. This was done in 
Krennrich ct al. (2001). By varying the cut tolerance 
t between 1.5 and 3.0 on a Crab Nebula dataset, the 
systematic error was determined to be ±0.05. 

The statistical error on the exponent of the flux con- 
stant is determined correspondingly VO.ll 2 + 0.10 2 = 
0.15 while the systematic error is derived by varying 
the absolute energy calibration by ± 15%. Note that 
the steepness of the spectrum makes this systematic 
error especially large. The final result is 

d£ _ 1n -7.3l±0.l5 st at±0.l6 syst 

x £-3.50±o.35 Bt at±o.05 By3t m _2 s _1 TeV 

3.7. Evidence for emission above 2.5 TeV 

The extended cuts are designed to minimize system- 
atic errors of the flux measurement. When purely sta- 
tistical evidence for the presence of emission is required, 
these cuts are inadequate for weak sources. This is the 
case for H1426 above 2.5 TeV, where the statistical sig- 




20 40 60 80 20 40 60 80 

ALPHA (°) ALPHA (°) 



Fig. 3. — As Figure 2, but only showing events with 
estimated energy above 2.5 TeV and using tighter cuts 
optimized on Crab Nebula data to give maximum sig- 
nificance above 2.5 TeV. The excess is expected at 
ALPHA< 12°. 

nificance of the excess obtained with the extended cuts 
drops well below 1 a in each energy bin. The question 
whether emission is present above a few TeV is, how- 
ever, so important for this source that a special effort 
is justified. 

Using 15.1 hours of contemporaneous Crab Nebula 
observations and a matching OFF dataset, we optimize 
a set of cuts to give maximum significance above 2.5 
TeV. With the standard extended cuts, the significance 
is 5.5ct. In the optimization, we vary independently the 
tolerance t for the upper cuts on WIDTH, LENGTH 
and ALPHA keeping the tolerance for the lower cuts on 
WIDTH and LENGTH constant at 2.0. In addition we 
require that at least 70 % of the excess obtained with 
the standard extended cuts is retained. The maximum 
significance (8.3 a) for the excess above an estimated 
energy of 2.5 TeV is reached for tolerance t = 0.5 for 
the upper WIDTH and LENGTH cuts and t = 1.9 for 
the ALPHA cut. These cuts are shown as dotted lines 
in Figure 1. 

When applied to the H1426 dataset (using a recalcu- 
lated energy estimation and correction K, see above), 
the optimized cuts improve the significance above 2.5 
TeV as one would expect in the case that emission is ac- 
tually present. Taking into account the error of the cor- 
rection factor K for the estimated energy of the back- 
ground events, there are 59 ± 22.7 excess events with 
estimated energy above 2.5 TeV. This corresponds to 



5 



(12±5) % of the excess rate observed from the Crab 
Nebula. The statistical evidence for emission from 
H1426 above 2.5 TeV is 2.6 a. Figure 3 shows the AL- 
PHA distribution of the H1426 events with estimated 
energy above 2.5 TeV after the optimized cuts. 

4. Discussion 

The spectral index found in this analysis agrees well 
with the first result (a = 3.55±0.5) published in Horan 
et al. (2002). The flux normalization constant is larger 
than that in Horan et al. (2002) but well within the 
statistical and systematic errors of this result which was 
obtained using a different analysis method and different 
Monte Carlo data. 

The only other published detection of H1426+428 at 
7-ray energies comes from HEGRA (Aharonian et al. 
2002) and is based on a dataset of similar size. The 
two results are compared in the following. 

From a total of 44.4 hours HEGRA obtains a to- 
tal number of excess events of 199.2 after loose cuts 
for their spectral analysis. The significance is 4.3 a. 
We obtain a slightly higher significance and an excess 
of 1540 events from 38.1 hours of data using similarly 
loose cuts. The relative numbers of events recorded 
by HEGRA and Whipple are consistent with the dif- 
ference in energy threshold which is a factor of w 2.5 
. An estimation of the integral spectral index from 
these numbers gives a value of 2.2 which already hints 
that the spectrum must be steep if there was no signif- 
icant change in the overall state of the source between 
1999/2000 and 2001. 

Direct comparison of the individual spectral points 
(Figure 4) shows a very good agreement indicating that 
there really was little change in the overall state of 
the source between the observing periods. Given the 
30% systematic errors of the absolute flux calibration 
of both our and the HEGRA measurement, this is con- 
sistent with the fact that the 7-ray rates measured by 
Horan et al. (2002) for the 2000 and the 2001 H1426 
dataset differ only by a factor 1.5. 

To make use of the total available information, we 
perform a fit of a power-law to the combined points 
from HEGRA and our analysis. No adjustments to 
the overall normalization of any of the datasets was 
performed since the points near 790 GeV agree very 
well. The errors on the Whipple points included in 
this fit take into account both the Gaussian error of 
the number of excess events (\/Aon + N Q g) and the un- 
certainty stemming from the energy estimation for the 




E (TeV) 



Fig. 4. — The differential energy spectrum of 
H1426+428 as measured in this analysis (filled points 
and upper limits) in 2001 and by HEGRA (open 
squares) in 1999/2000 (Aharonian et al. 2002). The 
84% confidence level upper limits were obtained using 
the method by Helene (1983). The solid line is the 
result of a power law fit to all points except the up- 
per limit at the lowest energy. The dashed line shows 
the scenario of the latter power law with an additional 
abrupt (super-exponential) cutoff still consistent with 
all data at the 95% confidence level. The minimum 
cutoff energy found for this scenario is 6 TeV. Also 
shown is the power-law fit to the HEGRA points only 
(dotted) and the power-law fit to the Whipple points 
only (dot-dashed). 
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background. The fit to all points (still excluding the 
upper limit at the lowest energies) results in a spec- 
trum which is essentially identical with what we obtain 
from our points only. However, the statistical errors 
are smaller: 

dF _ j^Q-r.seio.OTstat . £-3.54±o.27 stat m - 2 s - 1 TeV~ 1 
dE 

The reduced \ 2 01 this nt to H points is 0.94 . The 
measurements are consistent with the assumption that 
the spectrum of H1426 is a continuous powerlaw be- 
tween 250 GeV and 17 TeV. 

In order to further quantify the evidence for high- 
energy emission, we introduce a super-exponential cut- 
off in the measured power-law. This cutoff can be de- 
scribed analytically (exp(— (E/E ) 2 )) and is resonably 
abrupt such that it can serve as a parameter describing 
the energy above which there is no emission from the 
source. By keeping the nominal values of the power-law 
fit for normalization and spectral index and reducing 
the cutoff energy Eq until the x 2 of a fit to all signifi- 
cant points (i.e. excluding the points shown in Figure 4 
as upper limits) has increased to 14.1 (95% confidence 
level for 7 degrees of freedom) we obtain E — 5.5 TeV. 

Independently, we verify if this value for E is consis- 
tent with the fact that we observe from H1426 12±5 % 
of the Crab Nebula 7-rate above 2.5 TeV (see section 
3.7). Taking into account the larger spill-over effects 
caused by the finite energy resolution of our instrument 
and the steep spectrum of H1426 and furthermore the 
15 % uncertainty of our absolute energy calibration, we 
find that a 95 % confidence level lower limit on the in- 
tegral flux of H1426 above 2.5 TeV can be put at 3 % 
of the Crab Nebula flux. Using our measurement of 
the Crab Nebula flux with the same instrument, this 
corresponds to 

F m 426(E > 2.5TV) > 1.06 x lO^m-V 1 (95%c.l.) 

Varying the cutoff energy Eq introduced above until the 
integral flux above 2.5 TeV has the value of this lower 
limit, we find a value of Eq = 6 TeV which indepen- 
dently confirms the result from above. The spectrum 
with E = 6 TeV is shown as a dashed line in Figure 4. 

The HEGRA collaboration attempts to reconstruct 
the intrinsic energy spectrum of the source by assum- 
ing a model for the intergalactic IR background and 
calculating the expected absorption for 7-rays coming 
from H1426 (rcdshift 0.129). They arrive at an intrin- 
sic spectrum with index 1.9 and find that the "upturn" 
in the spectrum around a few TeV may be explained 



by a decrease of the slope of the 7-ray absorption as a 
function of energy. Our data is not inconsistent with 
the absorbed spectrum derived by HEGRA. In fact, 
our lowest point at 370 GeV agrees very well with the 
hypothesis of the absorbed spectrum, but not with the 
extrapolation of the alternatively fitted power law (dot- 
ted line in Figure 4) . The latter can be excluded with a 
confidence level > 99.5% which argues for a flattening 
of the spectrum at a few TeV. More data is needed to 
confirm this result. 

Between 250 GeV and 800 GeV, our data clearly 
favours a spectrum which is steeper than that of any 
other known TeV blazar at these energies. Since H1426 
is also the most distant known TeV blazar, the steep- 
ness of the spectrum may be interpreted as evidence 
for gamma-ray absorption in the intergalactic medium. 
However, the sample of TeV blazars still needs to be en- 
larged and the intrinsic spectrum well separated from 
absorption effects before one can come to firm conclu- 
sions. 
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